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Abstract 

We present a perturbative model for crystal-field calculations, which keeps into account the 
possible mixing of states labelled by different quantum number J. Analytical J-mixing results are 
obtained for a Hamiltonian of cubic symmetry and used to interpret published experimental data 
for actinide dioxides. A unified picture for all the considered compounds is proposed by taking 
into account the scaling properties of the crystal-field potential. 

PACS numbers: 75.10.Dg, 75.50.-y 
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I. INTRODUCTION 



Crystal-field (CF) theory^ is one of the most powerful theoretical methods to deal with 
the magnetic properties of rare-earth (RE) and actinide (An) ions, and Stevens' operator 
equivalents formalism is still the most commonly used to analyze experimental data due to 
its simplicity Unfortunately, this approach concentrates only on the CF splittings within the 
lowest-lying 2S+1 Lj multiplet of the considered ion, completely neglecting the contributions 
of excited multiplets ("J mixing"). Although the task of diagonalizing the large matrices 
related to the full f n configuration, including different J multiplets, is relatively easy to per- 
form numerically by means of today's computers, Stevens' approach often makes it possible 
to obtain analytical expressions for physical quantities of interest for systems of sufficiently 
high symmetry, thus leading to a deeper insight on the physics of several compounds. 

In the present paper, we discuss a perturbative approach which retains the validity of 
Stevens' formalism while correctly taking into account J-mixing effects. This method, which 
has led to interesting results for transition-metal (TM) based molecular clusters^ and fer- 
romagnetic exchange-driven RE-TM intermetallic compounds^ is now applied to evaluate 
the intramultiplet CF splittings in light An and RE ions. We exploit the method to an- 
alyze the CF of actinide dioxides. These large-gap semiconductors are among the most 
studied actinide compounds. Although /-electrons are well localized, the complexity of 
the magnetic Hamiltonian, which includes CF and magnetoelastic single-ion interactions, 
phonon-transmitted quadrupolar interactions, and multipolar superexchange couplings be- 
tween neighboring ions, leads to a number of interesting and unusual physical phenomena. 
Among them, we mention the proposed octupolar phase transition in Np02r*^ the observed 
CF-phonon bound states in Np02,— and the peculiar static and dynamic phenomena pro- 
duced by magnetoelastic interactions in U02 , 8l9ll ° some of which are not yet fully understood. 

The CF potential is the fundamental building block of any theoretical model of the prop- 
erties of dioxides, since this influences the single- ion behavior to a large extent. In particular, 
it determines which degrees of freedom of the / shell are left unquenched and the size of the 
corresponding multipole moments, which account for the low-T physical properties. Most 
of the published theoretical approaches are based on the above-mentioned Stevens' treat- 
ment of the CF, which includes only the lowest Russell-Saunders or Intermediate-Coupling 
multiplet of the ion. If one takes as starting point the CF of UO2, on which very detailed 
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information is available by inelastic neutron scattering (INS) experiments, 11 then scaling the 
CF of UO2 within the Stevens' framework (to take into account the different ionic radii) 
provides a good CF model for NpCv However, the same scheme applied to PUO2 is only 
qualitatively satisfactory, since it reproduces the correct level sequence but it underestimates 
the observed energy splitting. 

Moreover, this approach is not internally consistent because the so-obtained CF param- 
eters yield different results when additional ionic multiplets are included in the calculation. 
On the other hand, the increased complexity of J mixing calculations makes particularly 
hard to find CF parameter sets working consistently over the various compounds. Indeed, 
different sets have been proposed so far for dioxides within J mixing calculations. In par- 
ticular, in Np0 2 two distinct and equally good sets of parameters had been obtained. 12 

By our perturbative J mixing approach, we have been able to obtain a unique set which 
works well over all the considered compounds. 

II. THE PERTURBATIVE J-MIXING MODEL 

Following Ref. Q, the total free- ion and crystal- field Hamiltonian H = Hpi + Hcf-, with 



k.q 



can be rewritten in the form 



H = H + H 1 + H 2 , (2) 



where in the present case H coincides with Hpi and Hi and H 2 are chosen so that the 
former has nonzero matrix elements only between states belonging the the same 2S+1 Lj 
multiplet. It is possible to define a Hermitian operator Q such that the matrix element of 
the transformed Hamiltonian H' = e~ tQ He in are very small in the off-diagonal blocks, thus 
restoring the possibility to use an isolated-multiplet approach. In this framework,— 

(aJM \H'\ aJM') = E 0q j5 M m> + (aJM \Hi\aJM') (3) 
(aJM \H 2 \ a"J"M") (a"J"M" \H 2 \ aJM') 



En n ii jii — En n j 



IJII M I, 

where a, J label free ion manifolds and E 0a j are the eigenvalues of H . For clarity, in the 
following we will label the states as in the Russell- Saunders scheme, where a coincides with 



(L, S) and any additional quantum number necessary to identify the terms; yet the actual 
calculation will include Intermediate Coupling corrections to the eigenf unctions. 

Once we limit our calculations to the ground J multiplet only, the first term on the right- 
hand side of the above equation represents a uniform energy shift of the whole multiplet, 
while the second term is the usual ground-multiplet CF Hamiltonian H^ J \ The effect of J- 
mixing is accounted for by the third term, which will be considered as an extra contribution 
to the ground-multiplet Hamiltonian and labelled H^J X (we maintain the redundant super- 
script (J) notation in order to emphasize that the newly obtained J-mixing Hamiltonian 
also acts on the ground multiplet only). 

In the case of light actinides, the most important J mixing contribution comes from the 
two lowest J-multiplets, i.e. 2S+1 Lj and b+ L J+1 . Therefore, for the sake of simplicity we 
restrict our analysis to the case of these two multiplets, separated by an energy gap A by 
the spin-orbit interaction. From Eq. |3J H$fL can be written as 



JM 



H {J) 

mix 



M" k,q k',q' 



{JM 



(k) 



J + 1M") (J+1M" 



r (k') 



JM' 



A 

The Wigner-Eckart theorem, in the form 



(4) 



(J 1 M 1 |Cf >| J 2 M 2 ) = (-l) Jl ~ Ml (J x ||C (fc) || J 2 ) ( Jl kJ2 \ (5) 

\-M 1 q M 2 ) 

allows us to get rid of the sum over M" in Eq. HJ since the 3j symbol in Eq. El equals zero 
if M 2 7^ Mi — q. The products of 3j symbols can be rewritten as linear combinations of 
matrix elements of Stevens operators O q k ^ so that 

H^L= E (6) 

k,q,k',q' 

where the "mixing operators" M™ are^ 

n 

M™ = Y, c ti°v (7) 

p=0 



with conveniently defined c^> coefficients. 
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Let us consider the simple but important case of cubic symmetry, for which the CF 
Hamiltonian (0) has the form 



H, 



CF 



d 4) + 



a 



((>) 



(8) 



Restricting the calculations within the ground multiplet only and using the Stevens' operator 
equivalents formalism, Eq. (JBJ) becomes 



= ^ (01 + 50 4 ) + ^ 7 (Og - 2lOi) , 



B 



16 



(9) 



/3 and 7 being the fourth- and sixth-order Stevens factors. To introduce another common 
notation,— we define 

V. = f = A 4 <r 4 );y 6 = § = A 6 <r 6 >, (10) 

where (r n ) are the expectation values of the r n operator over the appropriate /-electron 
wavefunction. It is found that H^] x maintains the cubic symmetry, and has the form 

H%1 = (0\ + 50 4 4 ) + u 6 (Ol - 210 6 4 ) + u 8 (0° 8 + 280 8 4 + 650 8 8 ) + . . . , (11) 

where we did not explicitly write the terms containing operators of rank higher than 8 since 
01 = for k > 2 J, and J = 9/2 is the maximum possible value for the ground state of light 
lanthanides and actinides.— The coefficients appearing in Eq. (|11|) are dependent on J and 
can be written as 



A 



^ 4) (y 4 (J \\c^\\ J + i)) z + 4 6 ' 6) (v & (j \\c^\\ j + 1» 



(6,6) 



+^ 4 < 6) (y 4 (J \\c^\\ j+i)v,(j \\c^\\ j + 1»] , (12) 

and a list of v^ 1 ^ for the ground multiplets of f n configurations with 1 < n < 5 is given in 
Table [I] The reduced matrix elements (J ||C^|| J + l) are calculated by using the Inter- 
mediate Coupling free-ion wavefunctions. 

Ions with 6 / electrons have a J = ground singlet, so that no intramultiplet energy 
splitting can exist and J-mixing effects are evident only in the wavefunction composition; 
this case is then impossible to study by the present approach. As for ions with half-filled / 
shell such as Gd 3+ , Cm 3+ , and Bk 4+ , J mixing is generally negligible and it hardly affects 
any physical property. Finally, the perturbative J-mixing approach is in principle suitable to 



study heavy /-electron ions (8 < n < 13), taking into account that the ground J multiplet 
is mixed with J — 1 states, instead of J + 1. However, in this case the advantage of using 
the perturbative model with respect to the numerical diagonalization of H over a complete 
f n basis could be significantly reduced, because: i) the use of Stevens operators of rank 10 
(12) becomes necessary for J > 5 (6); ii) it is not always possible to diagonalize H^p + H^ 
analytically, even in cubic symmetry. In any case, the expected J-mixing strength is much 
smaller for heavy than for light elements. 



III. A PARTICULAR CASE: THE f 1 CONFIGURATION 

In order to show in detail how the present model can be applied to the study of rare-earth 
and actinide compounds, let us start from the simplest possible configuration: an ion with a 
single / electron. The only interaction present in the free-ion Hamiltonian is the spin-orbit 
coupling, 

Hpi = AL • S. (13) 

The f 1 spectra is composed of two multiplets only, 2 F^/2 (ground state) and 2 i 7 V/2; adding 
a crystal field of cubic symmetry, the 14 x 14 matrix representing the 2 F term is made up 
of two 4x4 and two 3x3 diagonal blocks. The complete Hamiltonian H can then be 
analytically diagonalized, and the resulting energy gap between the I\ doublet and a T 8 
quartet composing the ground multiplet is 

1 



Er 7 — E Ta 



1716 



3744\/ 4 - 4480\/ 6 



+ ^25600 (13^4 - 84\/ 6 ) 2 - 137280 (13V 4 - 84V 6 ) A + 9018009A 2 
-^16384 (13V 4 + 70V 6 ) 2 - 219648 (13\/ 4 + 70 V 6 ) A + 9018009A 2 



(14) 



Although the perturbative approach is not particularly useful here since an exact analytical 
solution can be obtained, let us study this simple case in order to clarify the details of the 
process and to understand its limits of validity. No intermediate coupling occurs, and (3 = 
2/315, 7 = 0, (J ||C 4 || J + 1) = 2V/KV77, and (J ||C 6 || J + 1) = -10^/2/143; moreover, the 
spin-orbit gap can be expressed as A = A( J+l) = 7A/2, with A = 100.5 meV. Diagonalizing 
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Hqp + H m l x (see Appendix for details), we find 

16 Tr If 20480 Tr2 614400 TrTr 204800 T .A , , 

Er, - E T , = —Vi + (— V? - — V lVe + ^Vi) , (15) 

which corresponds exactly to the series expansion of Eq. (JUjl up to the first order in A -1 . 

The intramultiplet energy gap calculated above can be directly measured by means of 
spectroscopic techniques; for example, inelastic neutron scattering measurements for PrC>2 
(where praseodymium ions have valence 4+, therefore presenting a 4/ 1 electronic configura- 
tion) have shown that E^ 7 — E^ 8 = 131 meV.— Figure Q shows the possible solutions of this 
equation in terms of the crystal- field parameters V4 and Vq, with three different expressions: 
the full black line corresponds to the exact diagonalization of the complete f 1 Hamiltonian 
[Eq. ([14p]: the dashed vertical line is obtained by neglecting J mixing and using Stevens' 
approximation E r? — E rs = (16/7)1^; finally, the dashed-dotted line takes J mixing into 
account perturbatively by the present approach [Eq. (fTKj) ]. The results are satisfactory: the 
agreement between the exact (full) and approximate (dashed-dotted) curves is qualitatively 
much better if compared with Stevens' approximation (dashed line), and the quantitative 
contribution of the excited states is correctly estimated in the small-CF range |T4/A| < 1. 
It may be noticed that the three curves of Fig. ^ almost coincide when Vq ~ 0, since in 
this case the largest part of the J-mixing correction (which is linear in Vq) vanishes [this 
can be verified by observing the relative magnitude of the different coefficients in Eq. (|15|)] ■ 
For the same reason, for small Vq, the value of V4 is slightly underestimated for Vq < 
and overestimated for Vq > 0. As the spin-orbit interaction is stronger for heavier ions and 
the gap between the two lowest multiplets grows with J, we expect our model to have a 
reasonably good performance over the whole lanthanide and actinide series. 

IV. THE CUBIC PHASE OF ACTINIDE DIOXIDES 

In this Section, the perturbative J-mixing model outlined so far will be applied to inter- 
pret the intramultiplet crystal-field splittings observed by inelastic neutron scattering (INS) 
for actinide dioxides AnC>2 (An = U, Np, Pu). The crystal- field analysis will be performed 
in terms of the parameters A 4 and Aq [Eq. (fTUjl ] instead of V4 and Vq, using the values of 
(r n ) given in Ref. 0. Although recent density functional studies have pointed out a certain 
degree of covalency for the An-0 bond in these systems,— this does not prevent one from 
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using the crystal field theory to analyze experimental data; it would be quite more difficult 
than in a ionic compound to calculate the CF parameters from first principles, but this is 
not our aim. We will demonstrate that the present method can be used for an analytical 
study of the experimental results over a wide range of parameters and compositions. 

INS spectra for UO2 in the paramagnetic phase (above = 30. 8K) display peaks 
between 150 ad 185 meV , 11 ' 18 which have been attributed to transitions between the r§ 
ground stat o 11 ' 18 ' 19 and excited 1^3 and I\ states, and no other magnetic transitions were 
reported up to 700 meV— (it is worth to recall that the T 5 — > Fx transition is not dipole- 
allowed so, if present, it will display a very small intensity with respect to the other two 
transitions). 

Figure 121 shows the values of A 4 and A$ for which the possible transitions lie within 
the experimentally observed range according to the perturbative model. In the previous 
Section, we have studied the quantitative discrepancy between the model's predictions and 
the exact results, which was found to be significant for large CF parameters. Following these 
estimates, we have determined a "safe zone" (indicated in Fig. El by a dashed ellipse) within 
which the true set of CF parameters for UO2 is located with high degree of confidence. 

For the paramagnetic phase of Np02 we have followed Amoretti et al.,— who observed 
a broad magnetic signal centered at 55 meV in the INS spectra and attributed this peak 
to a transition between the two r 8 quartets. They gave two possible solutions for the 
paramagnetic phase, labelled 3 and 4 in Fig. El we show that actually there are infinite 
possible solutions, divided in two branches. 

Pu02 displays a temperature-independend magnetic susceptibility below 1000 K,— so 
that the CF ground state is expected to be the T 1 singlet. Magnetic-dipole matrix elements 
involving this state within the multiplet are zero except with the triplet. Indeed, 
only one peak centered at 123 meV is observed in PUO2 INS spectra.— According to our 
perturbative model there are again infinite possible solutions, plotted in Fig. which cover 
a large area of the A^-vs-Aq diagram. 

It is clear that, in the case of An02, the INS data analysis cannot give unambiguous 
results if every compound is treated separately; this can also be inferred from the widely 
scattered sets of parameters which are found in the literature (some of which are listed in 
Table |n]and displayed in Fig. |3J). On the other hand, if we assume that the A^ parameters 
are approximately the same for all isostructural compounds,— an inspection of Fig. El shows 
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that the only area of the diagram where common solutions for An =U, Np, Pu might exist 
is around the point labelled 4, which corresponds to one of the two solutions (the "strong 
J-mixing" one) proposed by Amoretti et al.— for NpC>2 {A 4 = —19.6 meV/a^ ; A G = 
0.666 meV/ag). In order to verify this result, we have calculated the INS transition energies 
for UO2, NpC>2, PUO2, and PrC>2 with this set of parameters by numerical diagonalization 
of the complete f n configuration Hamiltonian ( Table llTT|) . Moreover, as a further test, we 
have calculated the magnetic susceptibility for Am02 and CfC>2 with the same parameters 



(Fig. the measurements can be found in Ref. |22j and Ref. |23j respectively). 

In all the cases examined so far, the comparison with experimental results is quite good, 
considering that a 100% exact scaling of the CF potential is not expected to hold. Therefore, 
our results point towards a coherent unified picture for the CF potential in actinide dioxides. 

The splitting predicted for PrC>2 seems less satisfactory as the measured value of the gap 
E(T 7 ) — E(T 8 ) = 131 meV is quite larger than the value calculated with the solution we 
propose. However, the Pr02 case is complicated by magnetoelastic interactions that are 
known to affect heavily the physics of this compound by increasing the bare value of the CF 
gap.— Hence, a value of this bare gap of the order of 80 meV is fully realistic. 

One feature which cannot be accounted for by the CF models proposed so far is the 
temperature-independence of the magnetic susceptibility of Pu02 above 600 K. Indeed, it 
is obvious that if there is a magnetic gap of 123 meV as observed by neutron scattering, 
when this level becomes thermally populated it will contribute to the susceptibility, no 
matter the mechanism which generates the splitting (unless this contribution is accidentally 
canceled by other contributions). Indeed, even a completely different approach based on 
Density-Functional-Theory calculations cannot remove the discrepancy between magnetic 
susceptibility and neutron scattering experiments.— Using the diagram in Fig. 01 as a guide 
we have been able to find a solution for Pu0 2 (A4 = —26.7 meV/ap, Vq = 1.68 meV/ag) 
which leads to Er 4 — Er 1 = 134 meV and a flat x(T) curve below 1000 K. For this set the 
Curie contribution of the excited triplet and the off-diagonal Van Vleck contribution of the 
Ti — T4 pair accidentally combine into an almost T-independent susceptibility. However, 
this solution is quite unstable even for small variation of the CF parameters; moreover, it 
does not give good results if applied to U0 2 and Np0 2 . 
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V. CONCLUSIONS 



We have developed a perturbative approach to J mixing which allows the Stevens' formal- 
ism to be recovered by replacing the original CF hamiltonian with an effective one operating 
within the same J multiplet. We have applied the method to the study of the CF in actinide 
dioxides. These compounds have extremely interesting physical properties, which are deter- 
mined by the CF to a large extent. It is therefore important to reach a precise understanding 
of the CF, which is the basic building block of all theoretical efforts devoted to understand 
these properties. Most studies of the CF potential rely on the Stevens' approach, which 
yields a solution working fairly consistently over several compounds. Yet, this solution does 
not work satisfactorily anymore when J mixing is included, and the practical impossibility 
to perform exact J mixing calculations in wide ranges of the parameter space has prevented 
so far the identification of a better alternative. By hugely decreasing the numerical effort 
in favor of analytical calculations, our method has allowed us to find a J mixing solution 
(A4 = —19.6 meV/ag , Aq = 0.666 meV/ag) which works consistently over all the considered 
compounds. 

Although in the present paper we have only studied cubic systems and included the 
two lowest multiplets only (which in dioxides account for almost 100% of J mixing), the 
method can be used for any point symmetry, at the price of additional terms in the effective 
Hamiltonian; also, as many multiplets as necessary can be included in the calculation, leading 
to additional contributions to the Uk coefficients of Eq. ()12j) . Even in this case our method, by 
allowing to perform a quantitative analysis of experimental data by means of a Stevens-like 
Hamiltonian containing higher-rank operators, is much more efficient than fully numerical 
diagonalization. In particular, wide ranges of the parameter space can be easily investigated, 
allowing to produce quite easily diagrams such as Figs. |2 and H3 which would otherwise be 
very hard (if not impossible, in some cases) to obtain. 
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APPENDIX: ANALYTICAL DIAGONALIZATION OF A CUBIC CRYSTAL- 
FIELD HAMILTON! AN 

This Appendix will be devoted to illustrate the results of analytical diagonalization of 
the cubic crystal-field Hamiltonian 

= j (Ol + 50t) + b f 6 (01 - 2101) + |L (0 8 ° + 280 8 4 + 650 8 8 ) , (A.l) 

which corresponds to Hq F + H^ x discussed in this paper [Eqs. (JHJ) and (|TT|) ] with 5 4 = 
8V 4 /3 + z/ 4 ; 6 6 = 16Ve7 + u 6 ; b 8 = z/ 8 . 

For J = 5/2 (Ce 3+ , Pr 4+ , Sm 3+ , Am 4+ , Pu 3+ . . .) 6 6 and b 8 have no effects, and the 
multiplet splits in a T 7 doublet and a r 8 quartet separated by E-p a — Er T = 456 4 . 

A J = 4 ground multiplet (Pr 3+ , U 4+ , Np 3+ , Pu 4+ . . .) is splitted into a singlet (Pi), a 
doublet (r 3 ) and two triplets (r 4 and r 5 ) by a cubic crystal field. The energy separations 
between these levels are: 

E Fa - E Tl = 180(-6 4 + 636 6 - 216 8 ); (A.2) 



E T4 - E Tl = 105(-6 4 + 636 6 - 2166 8 ); (A.3) 

E Tb - E Tl = 45(-96 4 + 1056 6 - 2806 8 . (A.4) 

For J = 9/2 (Nd 3+ , Np 4+ , U 3+ . . .) the multiplet splits into a T 7 doublet and two T 8 
quartets, respectively splitted by 

E ( 2) - £ ( x) = 15W7(103^ + 47046 4 6 6 + 856806§ - 105846 4 6 8 + 3628806 6 6 8 + 31298406§) 

8 8 V 

(A.5) 

and 

735 (-^r^ ~ E r (1) ) 

E F7 - E rW = —b A - 12600& 6 + 283506 8 + ^— a —. (A.6) 

1 8 2 2 

The above results, with their algebraical signs, are correct for any particular order of the 
energy levels. 
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TABLE I: («/) coefficients for values of J corresponding to the ground state of light lan- 

thanides and actinides. 



Coefficient 



J = 5/2 



J = 4 


J = 9/2 


— 1328/2760615 


— 296/920205 


-a/9/35 x 956/184041 


— a/17/7 x 198/184041 


—6464/920205 


—39488/1 5643485 


-1928/93648555 


-452/52026975 


-^2/35 x 32/212355 


4/ (212355VH9) 


-512/22297275 


-1072/126351225 


2/36891855 


4/184459275 


-a/2/35 x 2/250965 


-1/ (150579^119) 


7/418275 


49/14221350 



,(4,4) 
4 

,(6,6) 
4 

,(4,6) 
4 

.(4,4) 
6 

,(6,6) 
6 

,(4,6) 

,(4,4) 
8 

,(6,6) 



,(4.6) 



32/10395 
64/ (33VA55) 
448/6435 















TABLE II: Crystal-field parameters found in the literature for various actinide dioxides. The 
labels correspond to those used in Fig. |3J 



Label Compound Va (meV) V% (meV) Source 

1 U0 2 -409 24.8 Ref.ja 

2 U0 2 -123 26.5 Ref.Jl 

3 Np0 2 -104 6.2 Ref.L2 

4 Np0 2 -132 26.4 Ref.Jl 

5 Pu0 2 -151 31.0 Ref.21 
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TABLE III: The intramultiplet transition energies experimentally detected by INS for some ac- 
tinide and rare-earth dioxides are compared with the corresponding energy splitting calculated in a 
J-mixing framework with the unique set of parameters A4 = —19.6 meV/ag , Aq = 0.666 meV/ag 
(i.e. the solution labelled 4 in Fig. |3J). 



Compound 


Transition 


Exp. (meV) 


Energy (meV) 

Calc. (meV) 


uo 2 


r 5 - 




155 a 


167 


uo 2 


r 5 - 


+ r 4 


172 a 


187 


NpO-2 


r (2) 
1 8 


. r (i) 

"* 1 8 


55 6 


56 


Pu0 2 




+ r 4 


123 c 


112 


Pr0 2 


r 8 - 


+ r 7 


131 d 


82 



Q Refsrtfi! 

b Ref. . 
c Ref. n. 
d Ref. 15. 



«|andQ 
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FIG. 1: Possible solutions of the equation E-p 7 — E-p 8 = 131 meV in terms of the crystal-field 
parameters V4 and Vq of Pr0"2- Full line: J-mixing effects are taken into account by exact di- 
agonalization of the full f 1 Hamiltonian [Eq. (|14|) ]. Dashed line: J-mixing effects are neglected. 
Dashed-dotted line: J-mixing effects are taken into account perturbatively [Eq. (|15|) ]. 

FIG. 2: Possible solutions of the disequations: 165 meV < Ep 4 — Er 5 < 179 meV and 150 meV < 
Er 3 — Er 5 < 160 meV, in terms of the crystal-field parameters A4 and Aq of UO2. The dashed 
ellipse represents the "safe zone" defined in the text. 

FIG. 3: Possible solutions of the equations: -EL(i) — E ' ( 2 ) = 55 ± 5 meV (Np02); -Er 4 — Er x = 
123 ± 5 meV (PUO2), in terms of the crystal- field parameters A4 and Aq. The dashed ellipse near 
the bottom-right corner represents the "safe zone" for UO2 (see Fig. |^); the numbered dots are 
sets of parameters found in the literature for different compounds (see Table UTI for details). 

FIG. 4: Experimental (dots, Refs. and calculated (lines, present work) inverse magnetic 

susceptibility for Am02 and Cf02, with A4 = —19.6 meV/ag ; Aq = 0.666 meV/ag ; xXm (0) = 
80 mol/cm 3 ; XqI (0) = 3 mol/cm 3 . 
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